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ABSTRACT 


We  present  soft  7-ray  observations  by  the  Oriented  Scintillation  Spectrometer 
Experiment  (OSSE)  on  the  Compton  Gamma  Ray  Observatory  (GRO)  of  the 
transient  X-ray  binary  pulsar  A0535+26.  The  observations  were  made  1994 
February  8-17,  immediately  prior  to  the  peak  of  a  giant  outburst.  The  phase- 
averaged  spectrum  is  complex  and  cannot  be  described  by  a  single-component 
model.  We  find  that  structure  in  the  spectrum  above  100  keV  can  best  be 
modeled  by  an  absorption  feature  near  110  keV,  which  we  interpret  as  the 
signature  of  cyclotron  resonant  scattering.  Because  of  OSSE’s  45-keV  threshold, 
we  are  unable  to  make  a  definitive  statement  on  the  presence  of  a  55-keV 
absorption  line;  however  we  can  conclude  that  if  this  line  does  exist,  it  must 
have  a  smaller  optical  depth  than  the  line  at  110  keV.  A  first  harmonic  (= 
fundamental)  cyclotron  resonance  at  110  keV  corresponds  to  a  magnetic  field 
strength  at  the  surface  of  the  neutron  star  of  ~  1  X  1013  G  (~  5  X  1012  G  if  the 
first  harmonic  is  at  55  keV). 


Subject  headings:  gamma  rays:  observation  —  line  identification  —  stars: 
individual  (A0535+26)  —  stars:  magnetic  fields  —  stars:  neutron  —  X-rays: 


stars 
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1.  Introduction 

The  massive  X-ray  binary  A0535+26  is  a  recurrent  transient  pulsar,  with  pulse  period 
ss  104  s  (Rosenberg  et  al.  1975)  and  recurrence  period  ss  111  days  (Motch  et  al.  1991).  The 
optical  companion  is  the  Be  star  HDE  245770  (Li  et  al.  1979)  at  an  estimated  distance  of 
1.8  ±0.3  kpc  (Giangrande  et  al.  1980).  For  a  review,  see  Giovanelli  &  Graziati  (1992).  The 
X-ray  outbursts  are  thought  to  be  powered  by  enhanced  mass  transfer  from  the  companion 
star  near  periastron.  The  intensity  of  the  X-ray  emission  can  vary  dramatically  from  one 
periastron  passage  to  the  next:  No  outburst  may  be  observed;  a  “typical”  outburst  may 
occur,  with  peak  flux  of  a  few  hundred  mCrab  units  (1-10  keV);  or  a  “giant”  outburst  may 
occur,  with  peak  X-ray  flux  of  several  Crab  units.  Giant  outbursts  may  appear  away  from 
periastron  (Motch  et  al.  1991). 

Previous  hard  X-ray  spectra  of  A0535±26  have  been  modeled  variously  as  exponential 
( kT  =  17  keV;  Hameury  et  al.  1983),  blackbody  ( kT  =  9  keV;  Frontera  et  al.  1985), 
and  Wien  ( kT  =  8  keV;  Dal  Fiume  et  al.  1988).  Kendziorra  et  al.  (1992,  1994)  reported 
the  detection  of  two  absorption  features,  near  50  keV  and  100  keV,  in  the  spectrum  of 
A0535±26  in  outburst.  Previous  estimates  of  the  magnetic  field  strength  have  ranged 
over  4  —  20  X  1012  G,  based  on  the  suggestion  of  cyclotron  line  features  (Kendziorra  et  al. 
1992,  1994),  the  energy  dependence  of  the  pulsed  fraction  (Dal  Fiume  et  al.  1988),  or  the 
multiplicity  of  pulse  shapes  at  low  energies  (Kanno  1980). 

We  report  here  on  spectral  observations  of  A0535±26  by  the  Oriented  Scintillation 
Spectrometer  Experiment  (OSSE)  on  the  Compton  GRO  over  the  interval  1994  February 
8-17,  during  the  rise  immediately  preceding  the  peak  of  a  giant  outburst.  These  observations 
were  scheduled  as  a  “target  of  opportunity”  triggered  by  the  detection  of  the  outburst  by 
the  BATSE  instrument  on  GRO.  The  rise  to  maximum  of  this  giant  outburst  spanned  a 
binary  phase  interval  of  ~0.15,  with  the  peak  intensity  occurring  on  1994  February  18 
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at  binary  phase  0.95  ±  0.03,  according  to  the  ephemeris  of  Motch  et  al.  (1991).  At  the 
maximum,  the  pulsed  emission  reached  8  Crab  (nebula  +  pulsar)  flux  units  in  the  20-40 
keV  band  (Finger,  Wilson,  &  Hagedon  1994;  Wilson  et  al.  1994). 

2.  Observations 

The  OSSE  instrument  consists  of  four  nearly  identical  large-area  NaI(Tl)-CsI(Na) 
phoswich  detector  systems,  as  described  in  detail  in  Johnson  et  al.  (1993).  Full-resolution 
spectra  were  accumulated  every  4.096  s  in  a  sequence  of  two-minute  measurements  of  the 
source  field  alternated  with  two-minute,  offset-pointed  measurements  of  background.  The 
observation  was  made  in  twice-normal  gain  to  improve  the  low-energy  response.  A  total  of 
1.8  X  105  s  of  highest-quality  data,  equally  divided  between  source  and  background,  were 
collected  for  each  of  the  four  detectors. 

For  this  observation,  the  source  and  background  fields  were  placed  to  ensure  that  the 
Crab  Nebula,  a  potentially  contaminating  source,  did  not  cause  a  significant  contribution 
to  the  flux  measurement  of  A0535+26.  Below  100  keV,  the  net  (source  -  background)  Crab 
contribution  is  <0.1%  of  the  measured  flux  from  A0535+26,  and  in  the  critical  range  from 
100-200  keV  where  we  see  spectral  structure,  the  contribution  is  <  10-4  photons  cm-2  s_1 
MeV-1,  or  ^  1%  of  the  measured  A0535+26  flux. 


3.  Results 


We  detected  emission  from  A0535+26  with  high  statistical  significance  from  45 
keV  to  ~200  keV.  The  average  45-200  keV  luminosity  during  1994  February  8-17  was 
4  X  1036  (d/1.8  kpc)2  erg  s_1;  the  luminosity  doubled  over  the  nine  days  of  the  observation, 
from  2.5  x  1036  to  5.1  x  1036  (d/1.8  kpc)2  erg  s_1. 
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We  fit  the  phase-averaged  spectrum  between  45  keV  and  300  keV  from  the  sum  of 
the  four  detectors  for  the  entire  observation  using  a  forward-folding  method  (Johnson  et 
al.  1993).  We  derived  error  estimates  in  the  parameters  of  the  various  spectral  models  by 
%2-mapping,  with  confidence  intervals  specified  according  to  Lampton,  Margon,  &  Bowyer 
(1976).  With  such  a  strong  source,  it  is  important  to  address  possible  sources  of  systematic 
error.  We  have  identified  a  number  of  likely  contributors,  and  in  all  the  analysis  described 
herein  we  have  added  our  best  estimate  of  the  systematic  error  to  the  statistical  error.  The 
systematic  error  becomes  the  dominant  source  of  error  below  ss70  keV,  rising  to  ~3%  of  the 
flux  at  45  keV.  It  does  not  significantly  alter  the  values  of  the  best-fit  parameters. 

Fig.  la  shows  the  phase-averaged  spectrum  summed  over  the  entire  observation. 

The  spectrum  is  quite  soft  (Wien  temperature  sslO  keV),  although  no  single-component 
model  adequately  describes  the  spectrum.  Simple  continuum  models  show  very  significant 
residuals  and  have  unacceptably  large  %2  values  (Table  1).  The  structure  apparent  above 
100  keV  could  indicate  the  emergence  of  weak,  underlying,  high-energy  continuum  emission, 
or  it  could  indicate  a  broad  emission  or  absorption  line.  To  determine  which  of  these  is  most 
likely,  we  tested  a  number  of  models  in  two  classes:  a  class  of  two-component  continuum 
spectra,  and  a  class  of  single-component  continuum  spectra  together  with  line  profiles  for 
emission  or  absorption  (Table  1).  In  both  classes,  the  continuum  model  that  accounts  for 
the  bulk  of  the  emission  is  a  power  law  times  an  exponential  (hereafter,  “PLExp”),  chosen 
because  of  its  general  adequacy  in  describing  the  standard  and  hard  X-ray  spectra  of  X-ray 
binaries  (White,  Swank,  &  Holt  1983). 

The  models  from  the  first  class,  the  two-component  continuum  models,  do  not  give 
acceptable  fits  (Table  1).  The  weak,  hard  continuum  necessary  to  describe  the  emission 
above  100  keV — either  a  simple  power  law  or  PLExp — has  a  substantial  effect  on  the 
spectrum  below  100  keV  and  generates  significant  residuals.  Models  from  the  second  class 
employing  a  Gaussian  emission  feature  (centroid  energy  ssl50  keV)  are  similarly  poor,  with 
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%2  probabilities  <  10  7 . 

The  remaining  models  from  the  second  class,  those  with  broad  absorption  lines,  provide 
substantially  better  fits.  In  these  models,  the  continuum  is  modified  by  a  multiplicative 
transmission  function  from  one  or  more  absorption  lines,  which  we  presume  arise  from 
cyclotron  resonant  scattering.  First  proposed  by  Tanaka  (1986),  the  transmission  function 
assumes  an  exponentiated  Lorentzian-times-E12  form,  which  resembles  the  cyclotron 
scattering  cross  section  (Herold  1979;  Daugherty  &  Harding  1986;  Harding  &  Daugherty 
1991)  and  has  been  used  successfully  to  model  cyclotron  absorption  features  in  Ginga  data 
(e.g.  Nagase  et  al.  1991;  Tamura  et  al.  1992).  The  model  function  has  the  form 

f(E)  =  AC(E)eM~L(E)],  (1) 


where 


C(E)  =  E~a  exp (-E/kT) 


1  ’  h(E~  +  (nwy '  ( ’ 

The  continuum  term,  C(E),  is  the  PLExp  model  with  folding  energy  kT.  In  the  absorption 
profile  L(E),  parameters  r„,  n£x,  and  nW  are  the  optical  depth,  energy,  and  width  (the 
half-width  at  half  maximum  of  the  Lorentzian,  HWHM)  of  absorption  line  n  =  1  or  2.  We 
have  assumed  a  harmonic  relationship  of  1:2  between  the  line  energies  and  widths.  While 
in  the  relativistic  treatment  the  resonant  energies  are  not  strict  multiples  of  the  cyclotron 
energy  (Herold  1979),  the  nonlinearities  are  small  (of  order  several  percent  for  the  first  few 
levels)  and  have  been  ignored.  Trial  fits  with  independent  line  energies  and  widths  were 
consistent  with  the  harmonic  assumption  and  did  not  substantially  improve  the  fit. 


Table  1  shows  that  the  best  fit  is  achieved  with  the  general  PLExp  continuum  model 
with  two  absorption  lines,  although  the  first  harmonic  (=  fundamental)  near  55  keV  has  a 
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small  optical  depth  (ri  =  O.ltoii  j  68%  confidence,  Xmin  +  7.0),  which  serves  mostly  to  flatten 
the  model  spectrum  at  low  energies.  Indeed  this  absorption  line  is  not  strictly  required  with 
the  PLExp  continuum  model:  %2  probabilities  for  models  with  free  or  zero  optical  depth  at 
55  keV  are  both  acceptable  (33%  and  22%,  respectively).  The  95%  confidence  upper  limit 
is  r  <  1.1  (Xmin  +  12.6).  We  note  also  that  the  Wien  continuum  with  two  absorption  lines 
is  acceptable;  however  the  lines  are  extremely  broad  (42+3%  HWHM),  and  the  55  keV  line 
is  difficult  to  distinguish  from  a  smooth  change  in  the  continuum. 

Parameters  of  the  best-fit  model  are  given  in  Table  2.  The  errors  given  correspond 
to  68%  confidence  intervals  for  six  interesting  parameters  (we  judge  the  intensity  A  to  be 
irrelevant);  thus  we  use  Xmin  +  7.0.  Care  should  be  taken  in  interpreting  the  error  estimates 
for  this  non-linear  model;  for  example,  the  95%  (i.e.  “2cr”)  confidence  interval  is  not  given 
by  twice  the  68%  (i.e.  “ltr”)  confidence  interval.  The  best-fit  optical  depth  of  the  line  at 
110^4  keV  is  1.8^5,  implying  15%  transmission  at  the  centroid.  The  width  of  the  line 
(28+g6  keV  HWHM)  is  typical  of  the  absorption  features  detected  by  Ginga  (~  25%,  e.g. 
Tamura  et  al.  1992).  The  width  and  depth  of  the  line  combine  to  affect  the  overall  spectrum 
significantly  throughout  the  45-200  keV  band  (Fig.  lc).  Because  of  the  large  uncertainties 
in  the  data  above  150  keV,  the  existence  of  possible  higher  harmonics  near  165  keV  and 
220  keV  is  not  well  constrained;  e.g.  the  95%  confidence  upper  limit  on  the  optical  depth 
at  165  keV  is  r  <  4  (Xmin  +  14.1).  Fig.  la  shows  a  PLExp  continuum  model  with  a  single 
absorption  feature  (near  110  keV);  the  optical  depth  of  the  potential  first  harmonic  near  55 
keV  has  been  set  to  zero.  With  only  a  single  absorption  line,  significant  correlation  among 
the  parameters  is  removed,  and  the  confidence  intervals  (inset  Fig.  la)  are  significantly 
smaller  than  those  in  Table  2.  The  residuals  (Fig.  lb)  demonstrate  that  a  feature  near  55 
keV  is  not  required,  reinforcing  the  conclusion  drawn  from  the  %2  probabilities.  The  total 
excess  emission  above  the  model  in  the  200-600  keV  band  is  <2cr. 

We  have  also  employed  a  more  general  line  profile  than  that  of  eqs.  1  and  3,  a 
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multiplicative  Gaussian.  Fits  with  this  symmetric  profile  are  essentially  identical  in  quality 
to  those  given  above  (e.g.,  for  the  two-line  model  Xmin  =  1-09  with  89  degrees  of  freedom) 
and  lead  to  identical  conclusions,  in  particular  that  only  the  line  near  110  keV  is  strongly 
required.  The  line  centroid  is  shifted  to  a  slightly  higher  energy  (115+4  keV),  presumably 
because  the  Gaussian  lacks  the  blue  tail  of  the  exponentiated  Lorentzian-times-FJ2.  The 
line  width  is  unchanged  (26+r9  keV  HWHM),  although  it  is  not  well  constrained.  The 
line  width  is  governed  by  thermal  broadening,  as  well  as  the  averaging  over  magnetic  field 
conditions  and  viewing  angle  present  in  the  phase-averaged  spectrum. 

Fits  to  the  spectrum  for  each  day  show  that  the  line  centroid  and  width  vary  by 
<  2cr,  which  suggests  that  we  have  not  introduced  a  significant  bias  by  summing  the  data 
over  the  entire  observation.  However,  as  the  luminosity  increases  during  the  outburst,  the 
spectrum  does  soften.  In  order  to  characterize  the  spectral  change  in  a  simple  manner,  we 
fit  the  daily  spectra  over  45-80  keV  with  a  Wien  law.  Although  formally  these  fits  are  only 
marginally  acceptable  (%2  probabilities  typically  ~  1%),  they  are  sufficient  to  convey  that 
the  magnitude  of  the  softening  is  modest.  The  Wien  temperature  decreases  by  7%  as  the 
luminosity  doubles. 


4.  Discussion 

We  have  shown  that  an  absorption  feature  near  110  keV  is  present  in  the  phase- averaged 
spectrum  of  A0535+26  during  the  giant  outburst  of  1994  February.  This  is  the  first  time 
that  both  the  low-energy  and  high-energy  wings  of  an  absorption  line  have  been  clearly 
detected  at  such  high  energies  in  any  binary  pulsar.  We  interpret  this  absorption  as  a 
signature  of  cyclotron  resonance  scattering.  Two-continuum  and  continuum  +  emission-line 
models  are  strongly  rejected  by  the  %2  test  (Table  1). 

We  are  unable  to  make  a  definitive  statement  on  the  presence  of  a  55  keV  absorption 
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line  in  the  phase-averaged  spectrum.  Furthermore,  there  have  been  no  reports  to  date 
on  the  existence  of  such  a  feature  in  the  BATSE  data.  Adequate  fits  to  the  OSSE  data 
are  obtained  with  and  without  this  feature,  with  a  slight  preference  for  a  weak  line  in 
combination  with  the  general  power-law-times-exponential  continuum  (Table  1).  However, 
we  caution  that  it  is  always  difficult  to  distinguish  between  the  appearance  of  a  weak,  broad 
line  feature  and  an  inadequate  description  of  the  underlying  continuum:  There  is  significant 
coupling  between  the  line  and  continuum  model  parameters.  Such  is  the  case  here,  and  it 
is  compounded  by  uncertainties  in  the  instrument  response  and  the  absence  of  adequate 
spectral  coverage  below  the  nominal  line  energy.  We  can  conclude,  however,  that  if  a  55  keV 
feature  is  present,  its  optical  depth  must  be  significantly  smaller  than  that  of  the  110  keV 
feature.  It  is  surprising  that  the  ratio  of  depths  would  be  small  (~1:20  in  the  best-fit  model 
with  both  depths  free;  <1:3.5  at  95%  confidence).  Radiative  transfer  calculations  based  on 
a  neutron-star  paradigm  for  the  hard-spectrum  7-ray  bursts  do  show  significant  filling-in 
of  the  first  harmonic  by  two-7  decay  of  the  second  excited  state  (e.g.  Wang  et  al.  1989; 
Alexander  and  Meszaros  1991).  However,  Bulik  et  al.  (1992)  point  out  that  this  “photon 
spawning”  process  should  be  much  less  important  in  accreting  X-ray  pulsars,  where  the 
typically  soft,  thermal  spectrum  would  not  populate  the  higher  harmonics  significantly,  and 
therefore  would  not  lead  to  any  appreciable  downscattering  into  the  first  harmonic. 

Kendziorra  et  al.  (1992,  1994)  reported  the  detection  of  two  absorption  features  (ss50 
keV  and  sslOO  keV)  in  the  spectrum  of  A0535+26  at  periastron.  An  exponentially-truncated 
power  law  model  for  the  underlying  continuum  was  required  by  their  combined  TTM  and 
HEXE  spectrum;  a  Wien  law  was  unacceptable.  The  50-keV  absorption  was  detected  at 
low  significance  (^  3cr)  in  the  phase-averaged  spectrum,  and  HEXE  lacked  the  sensitivity  to 
detect  with  certainty  the  high-energy  wing  of  the  100-keV  line;  however,  most  importantly, 
the  amplitudes  of  the  absorption  features  depended  on  rotation  phase.  We  have  begun 
analysis  of  the  phase-resolved  spectroscopy  data  from  OSSE,  which  will  be  sensitive  to  a 
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phase-dependent  absorption  at  55  keV,  if  it  exists.  We  note  that  the  95%  confidence  limit 
we  set  here  on  the  optical  depth  at  55  keV  (r  <  1.1)  in  the  phase-averaged  spectrum  does 
not  exclude  the  HEXE  value. 

If  the  feature  at  110  keV  does  indeed  result  from  cyclotron  resonant  scattering  in  a 
strongly  magnetized  plasma,  it  can  be  used  to  estimate  the  magnetic  field  strength.  A 
first  harmonic  resonance  at  110  keV  corresponds  to  a  field  strength  at  the  surface  of  the 
neutron  star  of  9.5  X  1012  (1  +  z )  G,  where  z  is  the  gravitational  redshift.  If  instead  the 
first  harmonic  is  at  55  keV,  the  field  is  4.7  X  1012  (1  +  z )  G.  In  either  case,  this  is  the 
highest  field  strength  for  an  accreting  neutron  star  ever  measured  through  direct  detection 
of  cyclotron  absorption,  although  the  value  is  within  reason  and  is  consistent  with  the 
strongest  magnetic  fields  inferred  from  the  spin-down  of  radio  pulsars. 

This  work  was  supported  by  NASA  contract  S-10987-C  and  by  DARA  grant  50  OR 


92054. 
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Description 

Xdof  (dof) 

Single- continuum  models 

OTTBa 

59.4  (94) 

Wien 

11.7  (94) 

PLExpb 

5.44  (93) 

Two-continuum  models 

PLExp  +  Power  Law 

5.42  (91) 

PLExp  +  PLExp 

5.30  (90) 

Continuum  +  emission 

line  models 

Wien  +  Gaussian 

11.9  (91) 

PLExp  +  Gaussian 

4.72  (90) 

Continuum  with  absorption  models 

Wien  X  One  Linec 

1.86  (91) 

Wien  X  Two  Lines 

1.13  (90) 

PLExp  X  One  Linec 

1.11  (90) 

PLExp  x  Two  Lines 

1.06  (89) 

Table  1:  Minimum  Xdol  values  f°r  candidate  spectral  models. 

“Optically  thin  thermal  bremsstrahlung. 

bPower  law  times  exponential. 

cEl  ss  110  keV. 
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Continuum  Parameters  Line  Parameters 


Intensity  (cm'2  s_1  MeV-1)3 

o  oc+4.32b 
^.OJ_0.35 

Line  Energy,  1st  Harm.?  EL  (keV) 

55-Olli 

Power  Law  Index^  a 

—0.131°$ 

Line  HWHM,  1st  Harm.?W  (keV) 

Folding  Energy,  kT  (keV) 

17-8±|;| 

Optical  Depth,  1st  Harm.,  T\ 

O-lOloiio 

Optical  Depth,  2nd  Harm.,  r2 

1-81^ 

Table  2:  Best-fit  parameters  for  power-law-times-exponential  continuum  and  two  absorption 
lines. 

“At  70  keV. 

^Errors  correspond  to  68%  confidence,  assuming  six  interesting  parameters  (Xmin  +  7.0; 
Lampton,  Margon,  &  Bowyer  1976). 
ci.e.,  photon  spectral  index. 

^Energy  and  HWHM  of  2nd  harmonic  are  twice  those  of  1st  harmonic. 
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Fig.  1. —  a)  Phase- aver  aged  spectrum  of  A0535+26  from  1994  February  8-17.  The  outburst 
peaked  1994  February  18  (Wilson  et  al.  1994).  Solid  line  is  the  model  fit,  consisting  of  a 
power  law  times  an  exponential  continuum  (eq.  2)  with  transmission  function  of  a  single 
absorption  line  near  110  keV  (eq.  3).  Inset  shows  joint  68%  (%„;„  +  5.9)  and  95%  (%„;„  + 11.1) 
confidence  intervals  for  the  line  centroid  and  optical  depth  from  the  single-line  model,  b) 
Residuals  from  best  fit.  c)  Transmission  function  for  the  best-fit  absorption  line  near  110 
keV.  Data  are  divided  by  the  continuum  portion  (eq.  2)  of  the  best-fit  model. 


